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Abstract. In this paper, we first rework B. Kaufman's 1949 paper [8] by using representa- 
tion theory. Our approach leads to a simpler and more direct way of deriving the spectrum 
of the transfer matrix for the finite periodic Ising model. We then determine formulas for 
the spin correlation functions that depend on the matrix elements of the induced rotation 
associated with the spin operator in a basis of eigenvectors for the transfer matrix. The 
representation of the spin matrix elements is obtained by considering the spin operator as 
an intertwining map. We exhibit the "new" elements V+ and V- in the Bugrij-Lisovyy for- 
mula [2] as part of a holomorphic factorization of the periodic and anti-periodic summability 
kernels on the spectral curve associated with the induced rotation for the transfer matrix. 



1. Introduction 

In 1944, L. Onsager [12] determined the exact value of the specific heat as a function of 
temperature in the thermodynamic limit of the two-dimensional Ising model in the absence 
of an external magnetic field. He showed that the partition function can be approximated by 
the largest eigenvalue of the transfer matrix on the lattice. In 1949, B. Kaufman [8] simplified 
Onsager's calculation considerably by realizing the transfer matrix for the finite periodic Ising 
model as an element in a spin representation of the orthogonal group. The eigenvalues of the 
transfer matrix could then be found from the angles of rotation of the rotation matrix. In 
order to calculate the degree of order, the transfer matrix was diagonalized by transformations 
obtained as spin representations of orthogonal rotations. However, the calculation of the 
determinant of these rotations (the proof of whether they are proper or improper) is not 
included in Kaufman's paper. (See line 33, page 1237 in [8]). In this paper, we provide these 
calculations. 

There is a complex matrix T Z (V) whose entries are rational functions of z G C which 
completely determines the transfer matrix V on the finite periodic lattice. The matrix T Z {V) 
is called the induced rotation associated with the transfer matrix. The set of pairs (A, z) 
such that det(A — T Z (V)) = is an elliptic curve M. which is important for the spectral 
analysis of the transfer matrix. In particular, the map Ai 3 (A, z) — > z G P 1 is a two fold 
covering, and there are two cycles Ai± on Ai which cover the circle S 1 = {z : \z\ = 1}. On 
the cycle M.+ we have A < 1, and on the cycle we have A > 1. Just which points 

Zj G S 1 are relevant for the spectral analysis depend on the boundary conditions for the 
model. For spin periodic boundary conditions on the lattice, the (2M + l) th roots of unity, 
z 2M+1 = 1, are relevant as are the (2M + l) th roots of —1, z 2M+1 = —1. We will refer to 
these two finite sets as the periodic spectrum Ep and the anti-periodic spectrum S^. In the 
infinite-volume limit all the points z G S 1 are relevant. Kaufman [8] showed that the space 
in which the transfer matrix acts, can be divided into two invariant subspaces, which we will 
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denote by (U = 1) and (U = — 1). Here U is a product of the basis elements of the finite 
sequence space, W := l 2 {— M, M, C 2 ), which are certain representations of the Clifford 
relations. The expression (U = ±1) is the short-hand notation for the ±1 eigenspaces of U. 
More specifically, Kaufman showed that the transfer matrix V can be written as the direct 
sum V = V A © V p , where V A = V\( U=1 ^ and V p = V\( U= _iy Here we have introduced 
the letters A and P to refer to the restriction to the anti-periodic and periodic spectrum 
respectively. We rework Kaufman's paper [8] by proving that (U = 1) is unitarily equivalent 
to the even tensor algebra of a subspace of W in the anti-periodic Fourier representation 
while (U = —1) is unitarily equivalent to the even tensor algebra of a subspace of W in the 
periodic Fourier representation. This result is summarized in Theorem 3.1 and is the main 
result of this paper. In the anti-periodic Fock representation, we can choose a representation 
of the transfer matrix V A such that the vacuum vector Oa is an eigenvector associated with its 
largest eigenvalue. In a similar fashion, we can choose a representation of the transfer matrix 
V p in the periodic Fock representation such that the vacuum vector Op is an eigenvector 
associated with its largest eigenvalue. By using these representations of the transfer matrix, 
we can easily compute its full spectrum. We believe that this approach is much more direct 
and simpler than the method employed in Kaufman's paper. By applying this representation 
theoretic approach, we can also determine formulas for the spin matrix elements, in a basis 
of eigenvectors for the transfer matrix, that depend on the matrix elements of the induced 
rotation associated with the spin operator. The spin matrix elements will be calculated by 
regarding the spin operator as an intertwining map for the periodic and anti-periodic Fock 
representations. 

The paper is organized as follows. In Section 2, we describe the transfer matrix approach 
to the Ising model as given in [8], and we provide the formula for the induced rotation 
associated with the transfer matrix. In Section 3, we prove Theorem 3.1 which we described 
above, and apply it to derive the spectrum of the transfer matrix in Theorem 3.3. In Section 
4, we compute the matrix elements of the induced rotation associated with the spin operator 
by restricting the spin operator to be a map from the Fourier space with periodic boundary 
conditions to the Fourier space with anti-periodic boundary conditions. This result is given 
in Theorem 4.1. Let us denote the matrix of the induced rotation for the spin operator by 

s-=( A B ) 
s ' V C D ) ' 

where A, B, C, D are matrix elements for s in a polarization of W. The correlation func- 
tions on the cylinder and the torus can be evaluated in terms of spin matrix elements in 
an orthonormal basis of eigenvectors for the transfer matrix. We recall this calculation in 
Section 5. Subsequently, in Section 6, Theorem 6.1, we provide Pfaffian formulas for the spin 
correlation functions that depend on the inverse, D^ 1 . To invert D is an extremely diffi- 
cult task, and we only have a formula for it in terms of a Bugrij-Lisovyy conjecture for the 
spin matrix elements. Bugrij and Lisovyy proposed the explicit formulas for the spin matrix 
elements on the finite periodic lattice for the isotropic case [2] and for the anisotropic case 
[3]. We include these formulas in Section 7. A proof of these formulas was given in [5] and 
[6] but it is rather complicated. In Section 8, we express the proposed formulas for the spin 
matrix elements in terms of a product of Jacobian elliptic functions in the uniformization 
parameter of the spectral curve. In Section 9, we discuss the numerical comparison of our 
D~ l matrix elements with the corresponding term in the Bugrij-Lisovyy formula. In Section 
10.7, we turn to a calculation that shows the role played by the "new "elements V± in the 
Bugrij-Lisovyy formula in a factorization of the ratio of summability kernels on the spectral 
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curve. We summarize this result in Theorem 10.1. In Appendix A, we introduce the Fock 
representation of the Clifford Algebra of W as given in [13]. In Appendix B, we introduce 
a Berezin integral representation for the matrix elements of the Fock representation of an 
element in the Clifford group. We express the kernel of this operator as an exponential of a 
skew symmetric matrix whose matrix elements depend on the matrix elements of its induced 
rotation. We give this result in Lemma B.l and Theorem B.2. Finally, in Appendix C, we 
provide the formula for the spectral curve associated with the induced rotation for the trans- 
fer matrix. We recall some properties of the Jacobian Elliptic functions which are involved 
in the holomorphic factorization of the summability kernels on the spectral curve and in the 
product formula for the spin matrix elements. 

2. Transfer Matrix 

We begin this section by recounting the analysis in the 1949 paper of Bruria Kaufman 
[8] which shows that the transfer matrix for the finite periodic Ising model can be expressed 
as an element in a spin representation of the orthogonal group. We then continue with a 
calculation of the induced rotation associated with the transfer matrix. 

For positive integers, M and N, consider the set Im = {— M, — M + 1, M} and the finite 
lattice A := Im x In- Let each vertex be assigned a spin value of +1 (spin up) or —1 (spin 
down). A configuration a is a particular assignment of spin values to the vertices, i.e. a spin 
configuration is a map, 

a : A^ {+1,-1}. 

In this paper, a satisfies the periodic boundary conditions a(—M,j) = a{M + and 
a(j,—N) = cr(j,N + 1) for all j. Each spin interacts ferromagnetically with its nearest 
neighbors; in a configuration a the interaction energy in the absence of an external magnetic 
field is defined by 

E\(a) = — Jij&i&j- 
{i,j)e\,\i-j\=i 

The sum is over all nearest-neighbors in A with real valued interaction constants, Jij = 
J\ > 0, if the sites are horizontally separated, and Jy = J2 > 0, if the sites are vertically 
separated. The probability of a given configuration a is proportional to the Boltzmann 
weight, 

w(a) := exp ( - , 

where T is the temperature and ks is the Boltzmann constant. The total weight is given in 
terms of the partition function, 

Z A = ^ w(a), 

where S7a is the set of all possible configurations on A. The correlation functions are the 
expected values of products of spin variables at sites ji, ■■■,j n in A, 

• • • • CJ Jn>A = 7^" °J1 " ' • CJ in U; ( CT )- 

Let fiA(row) denote the space of configurations of a row. An i th row configuration is a map, 



< 7 i :/ Jltf ->{+l,-l}. 
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For i = —N, ...,N, we denote a collection of configurations a 1 G ^^(row) as cr*- := a^y Thus, 
the spin variable <Jij is located at the site j in the i th row. For a, r G J7a(tow) we define the 
2 2M+1 dimensional matrices, 

M M 

V\{o) := exp( ^ /CiCTjCTj+i), V 2 (o",t) := exp( ^ K. 2 OjTj), (2.1) 

j=-M j=-M 

where /Q := for Z = 1,2. Kaufman [8] showed that the partition function can be written 
as the trace of the 2N + 1 power of the transfer matrix V\ V2, 

Z A = tr[(V 1 V 2 ) 2N+1 }. (2.2) 

Introduce the matrices, 

1 \ /o 1\ / 1 

N -1 ) ' V 10 / V 01 

and define 

Oj := J • • • J 0cr J • • • I (2.3) 

M+j 

Cj := 7 • ■ ■ /0C0 /0 • • • 0/ 

M+j 

which act on the tensor product, 

M 

C 2 := C 2 M • • • C 2 M = C 2 ( 2M+1 ), 

j=-M 

where C 2 = C 2 for each j. One can pick a basis for this tensor product space such that the 
action of the spin operator <jj for a G ^(row) is given by (2.3) [13]. For — M < k < M, 
define 

£> fe := C • • • C 0o- / • • • /, 
y v ' 

M+k 

Qk ■= C • • • C - iaC 1 • • • I 

M+k 

which satisfy the generator relations for the Clifford algebra, 

PkPl + PlPk = 2(5 fe i, q k qi + qiq k = 25 M , p k q t + q x p k = 0. (2.4) 

Introduce the finite sequence space, W := Z 2 (/m;C 2 ), with an orthonormal basis 
^^k=-M sucri ^at for an element i> G W, we have 

M 

u = 75 x k(. v )qk + yk(v)p k . 

k=-M 

We define the distinguished nondegenerate complex bilinear form (-, •) on W by 

M 

(u,t>) = ^ x k {u)x k {v) + y k {u)y k {v). (2.5) 

fc=-M 



The Hermitian inner product on W is given by (u,v) = (u,v), where the conjugation is 
defined as 

M 
k=-M 

The representations {qk,Pk} are generators of an irreducible ^representation of the Clifford 
algebra, Cliff (W), on C 2 ( 2M+1 ) [13]. Now define the operator, 

M M 



U := ] | C k = Yl iPkW,- 

k=-M k=-M 



which plays a central role in Kaufman's analysis. It is shown in [8] that the transfer matrices 
V\ and V2 in (2.1) can be written 

v i = (lljl-M ex P (-ifciPj+iQj)) ex P (HCip- M qMU), 

Vl = UfL-M ex P ( i!C *2Pj<lj), 

where the dual interaction constant /C| is defined by the relation 

sinh(2/C2)sinh(2/C 2 ) = 1. 

It is easy to check that U 2 = 1 and that U commutes with even elements in the Clifford 
algebra. From this it follows immediately that V\ and V2 leave invariant the ±1 eigenspaces 
for U. Let (U = ±1) denote the ±1 eigenspaces for U. Define 

V := VxV2 

and 

M \ / M-l 



V A : = i(7 + U) ( Y[ expiK,2PjQj) ( Yi exp ) ex P (^iP-mQm), (2.6) 

, M \ f M \ 

V P : = -(I-U)( [] exp(HC* 2Pjqj ))( ]J exp (-HC lPj+iqj ) J . (2.7) 

The letters P and A refer to periodic and anti-periodic respectively, and the reason for their 
appearance will be clear shortly. A central result in [8] is that the transfer matrix V can be 
written as the direct sum, 



where 



V = V A © v p , 



V A = V\ {U=1) and V p = V\ (u= _ l) . 



The exponential factors in V A and V p are elements in a spin representation of the orthogonal 
group [8]. 

A linear transformation V on (g>jL_ M C 2 is an element of the Clifford group if there exists 
a linear transformation, 

T(V) :W -^W, 
such that for all w 6 W C Cliff (W), we have 

VwV- 1 = T{V)w. 
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The operator T(V) is called the induced rotation associated with V. The induced rotation 
T{V) is complex orthogonal with respect to the bilinear form (•, •) and it determines V up 
to a scalar multiple [13]. For j = 1, 2 introduce the short-hand notation 

cj := cosh(2/Cj) and Sj := sinh(2/Cj) 

for the hyperbolic parametrization of the Boltzmann weights. We write 

c* := cosh(2/C*) = ^ and s* := sinh(2/C*) = — 

for the dual interactions. Adapting the technique introduced in [13], we compute the induced 
rotation associated with the transfer matrix in (x, y) coordinates, 



T(V 2 



X \ = / c\ -%a\ \ / X 

y \ is *2 c 2 J V v 



X \ I C\ IS\Z \ I X 

y ) ~ v -^i^ c i / v y 



The second equation is understood as follows. The operator z acts on a finite sequence x 
by sending x& to Xk-i- On the invariant subspace for V\, where {U = —1), the sequence 
is extended to be 2M + 1 periodic. On the invariant subspace for V\, where (U = 1), the 
sequence is extended to be 2M + 1 anti-periodic. For k = —M, M, introduce 

3 P._ ^__27T(fc + i) 



and 



fe ' 2M + 1' fc ' 2M + 1 



z P := z P (k) = e ld * , := ^(fe) = e* e *. 



Define the two sets, S P = {z G C : z 2M+1 = 1} and E A = {z e C : z 2M+1 = -1}, which 
we refer to as the periodic and anti-periodic spectrum respectively. By specializing the finite 
Fourier series, 

1 M 

to the periodic spectral points, z £ Sp, on the subspace (U = —1) or to the anti-periodic 
spectral points, z £ T,a, on the subspace (U = 1), we obtain that T(V\) acts as a multiplica- 
tion operator, 

T(V\ ( X ^ \ - ( Cl iSlz l \ ( X ^ 
[Vl) { y(z) J ~ { -is lZ d J V V(z) 

Define the induced rotation for the symmetrized transfer matrix, 

T(V) :=T(y/)T(lq)T(y/). 
Then for temperatures below the critical temperature Tq we find, 

T(V) = ( _ cosh ^( z ) w(z)smhj(z) \ ^ g . 

y w(z)sinhj(z) cosh 7(2) J ' 

where 7(2) is defined as the positive root of 

cosh 7(2) = c^ci - s*si (^±§^) , (2.9) 
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and where 

Here Aj(z) is denned by 
where 



.Alc*- 1 )^*- 1 )' 

Aj{z) := - i for j = l, 2, 



ai := {cX-s\){c 2 + s 2 ) = e 2 ^- K *\ 
a 2 := ( c * +s *)( C2 + S2 ) = e 2 (^+^9. 
For T < Tc, we have that 1 < ai < a 2 , and the square root is chosen to lie in the right 



half plane with positive real part. We observe that z — > {Aj{z)) ± is analytic for z in a 



neighborhood of the unit disc {\z\ < a±) while z — > (Aj(z is analytic in a neighborhood 
of the exterior of the unit disc (\z\ > a^ 1 ). 

3. The Spectrum of the Transfer Matrix 

Kaufman [8] showed that the eigenvalues of the transfer matrix can be found from the 
angles of rotation of the rotation matrix by realizing the 2 2M+1 dimensional transfer matrix 
as a spin representation of the (2M + 1) dimensional rotation matrix. In this section we 
compute the eigenvalues for the transfer matrices V A and V p by using representation theory. 
This approach is simpler and more direct than Kaufman's method. 

After Fourier transform, the Hermitian inner product on W becomes 

M 

^2 x(z k )x(z k ) + y{z k )y(z k ) 

k=-M 

and the complex bilinear form becomes 

M 

x{z k )x'{ Z ^) + y{z k )y'{z k l ), (3.1) 

k=-M 

where the conjugation is given by x(z) i-> x{z~ 1 ). Let T A denote the induced rotation 
associated with V A , where V A is given in (2.6). Let T p denote the induced rotation associated 
with V p , where V p is given in (2.7). Both T A and T p have positive real spectrum which 
does not contain 1. We are interested in the isotropic splittings, 

W = W A and W = W^@ W p , 

where W A and W p are the span of all eigenvectors of T A and T p respectively associated 
with the eigenvalues between and 1, and W A and W p are the span of all eigenvectors of 
T A and T p respectively associated with the eigenvalues greater than 1. Define 



a(z) := 



/ A 1 (z)A 2 (z) 
Aiiz-^Aiiz- 1 ) 



which is normalized such that a(l) > 0. The eigenvalues of T p are given by e _7 ^ Zp ^^ 
and e 7(Ap ( fe )) for k = —M,...,M, where 7(2) is defined as the positive root of (2.9). The 
corresponding eigenvectors are 

= T2 ( -tr 1 ) s ^> {z) (3 ' 2) 
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and 



» == ^ ( _ jj^-i ) (3-3) 



1 if z = zp(k); 
ifz^zp(jfe). 



respectively for z G Sp. Here 

<** P (fc)(*0 = 

Then {e^ fc , fc }^£_ M is an orthonormal basis for W with respect to the Hermitian inner 
product. We also have (e^ fc ,e_j) = tffcj, so {e^ fc }^£_ M and {e^ k } 1 jf = _ M are dual basis for 
Wj^ and W p with respect to the complex bilinear form. The eigenvalues and eigenvectors, 
e A k (z), of T A are defined in the exact same way for z G S^. From (2.8) we see that VF_j5 
are the ±1 eigenspaces of the polarization given by the multiplication operators in the anti- 
periodic Fourier representation, 

«-) = -(w(L, M 

and Wl are the ±1 eigenspaces of the polarization given by the multiplication operators in 
the periodic Fourier representation, 

«">-(*(?,) 1 P) )' < 3 - 5 > 

Let Q A denote multiplication by Q(za), and let Q p denote multiplication by Q(zp). Define 

Q A :=\{I±Q A ) and := \{I ± Q p ) . 

Then 

w£ = Q A W and = Q P W 

so and Q± are orthogonal projections on W± and W± respectively. Define the alternating 
tensor algebras, 

Alt even (W + ) := Alt 2k (W + ) and Alt odd (^+) := Alt 2k+1 (W + ), 

0<2fc<n 0<2fc+l<n 

where Alt fc (W + ) is the space of alternating k tensors over W+, Alt°(W+) = C and n = 
dim(W+). Let F p and F A denote the Fock representations associated with the Clifford 
relations acting on Alt(W+) and Alt(W+) respectively. We have for x € W, 

F p (x) = c(x+) + a{x-), (3.6) 

where x± = Q±x, and c(-) and a(-) are creation and annihilation operators as defined in 
(A.4) and (A.5). Recall that V = V A © V p , where V A = V\ {u=1) and V p = %/=-i), and 

M 

where U = | J ipkQk- We will prove the following theorem. 

k=-M 

Theorem 3.1. Consider the isotropic splittings W = W A © W A and W = W p © W p 
associated with the polarizations defined in (3.4) and (3.5). Let (U = ±1) denote the ±1 

M 



eigenspaces for U , where U = j j ipkQk- Then 

k=~M 

Alt eve „(W^) ~ (U = 1), Altev«n(Wf ) ~ (17 = -1), 



Alt odd (W^) ~ (U = -1), AltoddCWf ) ~ (U = 1), 
where ~ denotes unitarily equivalent. 

Proof. The strategy of the proof is as follows. We first make a change of basis {^|, ^}^1_ M 

to a basis { -^^(^(A;)), t^p(- z (^))}^_j^-! where q(z(k)) and p(z(k)) are real with respect to 
the conjugation, v(z) i-> for z in Sp or in S^. The elements q(z(k)) and p{z{k)) will 

be defined below. Then define the "volume elements" , 

M M 
U A := H ip(z A (k))q(z A (k)) and U P := \{ i P (z P (k))q(z P (k)), 
k=~M k=~M 

in the Clifford algebra. Let M p and M A denote the number operators in Alt(W+) and 
Alt(Wjj5) respectively, which are defined as 

M p \A\t k {W p ) = k and J\f A \ A\t k (W A ) = k. (3.7) 

We show below that 

U P = {-1) NP and U A = (-l) NA . (3.8) 

Denote the linear transformation that sends ^= to - l =q(zp(k)) and ^= to ^p{zp(k)) by R p , 

and the linear transformation that sends ^ to ^q(zA(k)) and ^ to ■ 1 ~p( K z A (k)) by i?" 4 . 

The transformations, R p and -R" 4 , are real and orthogonal, and it follows from page 153 of 
[15] that they induce an automorphism of the Clifford algebra such that 

U = det(R p )U p and U = det(R A )U A . (3.9) 

We show below that 

deti? p = -l and deti? A = l. (3.10) 
Then it follows from (3.8), (3.9) and (3.10) that 

U = -{-l)* P on A\t(W p ) (3.11) 

and 

U = {-l) NA on A\t(W A ). (3.12) 

Since (-1)^ = 1 on Alt cvcn (W P ) and {-1) Ma = 1 on Alt cvc „(W^), we have from (3.11) and 
(3.12) that 

AltevBn(Wf ) ~ (U = -1) and Alt even (W A ) ~ (U = 1). 
Subsequently, it follows that 

Altodd(Wf) ~ (U = 1) and A\t odd (W A ) ~ (U = -1). 

We now return to the introduction of the basis elements, {-^q(zp(k)), -^p{zp(k))}fL_ M and 

then proceed to prove the first identity in (3.8). In the polarization, W = W p © W p , define 
the creation operators as 

p* ( p \ 

and the annihilation operators as 

P i P \ 



where c(-) and a(-) are defined in (A. 4) and (A. 5), and {e£ k } are the eigenvectors for the 
induced rotation T p as defined in (3.2) and (3.3). In the zp coordinates we define the basis 
elements for W as 

p(z P (k)):=(a P *+a P ) and q(z P (k)) := i(a P - a P *). (3.13) 

By using the anti-commutative relations given in (A. 6), it can be checked that p(zp(k)) and 
q{zp{k)) satisfy the Clifford relations. Since e p ^{z" 1 ) = e p k (z), it is not hard to see that 
p(zp(k)) and q(zp(k)) are real with respect to the conjugation, v(z) i-> v(z~ l ). It can be 
checked that (U p ) 2 = 1 and 

U P p(zp(k)){U P y l = -p(z P (k)) and U P q{z P {k))(U P )- 1 = -q(z P (k)), 

so U p is an element of the Clifford group with induced rotation —1 on W p . Furthermore, 
from (3.7) we have 

{-l) MT F p {v){-l) MV = -F p {v) for v e A\t(W P ). 

It follows that F P {U P ) and (-1)^ have the same induced rotations in the Fock representa- 
tion. Since we also have (-l) 2AfP = 1, it follows that F P {U P ) = h{-l) NP for h = ±1. We 
show that l\ = 1. From (3.13) we have that 

a k* = hiP( z P( k )) + Mzpik))) and of = \{p{z P {k)) - iq{z P (k))). 
In the Fock representation associated with the polarization above, the number operator N v 

M 

is given by M v = ^ a P *a P . Using the fact that (a P *a P ) 2 = a P *a P , following a similar 

k=-M 

argument given in [13], we obtain 

M M M 

{-l) NV = J] e™<*<= n (l-2a P *a P )= ]J ip{z P (k))q(z P (k)) = U P . 
k=-M k=-M k=-M 

The basis elements {-^q(zA(k)), -^p(zA(k))} in the za coordinates are defined in an analo- 
gous way, and the second identity in (3.8) can be proved in a similar fashion. We now return 
to the proof of (3.10). We start by proving that det R p = —1. The linear transformation R p 
consists of a compositions of transformations, 

R P . £2(2M+^ pS>Jrp > C 2(2M+1) Rl > C 2(2M+1) R ' 2 > ^2{2M+1) X — ^ q2(2M+1) . 

Here Tp is the finite inverse Fourier transform in the periodic representation, R p is the 
transformation from the Fourier series representation to {e p k , e p _ k } k v L_ M , while R P is the 
transformation from {a k *, a p k } k v L_ M to {q{zp(k)),p(zp(k))} k L_ M . 

Since {q(zp(k)),p{zp(k))} is not the normalized basis for W, but | -^=q(zp(k)), -^p{zp{k)) 

is, the transformation x i— >■ \/2x takes the first basis to the second. Let us denote the matrix 
of this transformation by R3. (Note that the factor of 2 is not incorporated into our definition 
of the Clifford relations: xy + yx = (x, y)e for x, y G W, where the more standard definition 
is xy + yx = 2(x, y)e.) It follows that 

det( J R 3 ) = 2 2M+1 . (3.14) 

The inverse Fourier transform 

JT p . C 2Af+l ^ C 2M+1 
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in the zp representation has the following matrix, 



T P = 



2nikl 



1 



V2M + 1 



-M 
-P,-M 



.. .. Z 



-M 
P,M 



i M M I 

\ Z P,-M ■■ ■■ Z P,M I 



where z l Pk := z l P (k) = e 2M + 1 for k, I = —M, M. It follows from [11] that the determinant 
of Tp is given as 



det(J» 



\M 



(3.15) 



/ x- M \ 

v-m 



Let Tp © Tp act on the vector 



It can be checked that by interchanging the rows 



XM 
\VM J 

and columns an even number of times in the matrix of Tp © Tp that 



det(Tp(BT P ) = (-1) M 



(3.16) 



Define for k = -M, ..,M, 



J_ ( a(z P {k)) 



a(z P (k)) 



Rl ' k ' ^2 V izp{k)a(z P {k))' 1 -iz P {k)a{z P {k))- 1 



The matrix of Rf is given by 



M 



k=-M 



where 



detR? =((-z) 2M+1 J] 

\ 7 A/T / 



■2M+1 



(3.17) 



k=-M 

The matrix of i?f that goes from {a^,a-k}%L_ M to {q(zp(k),p(zp(k)}^L_ M is given by 



( -i 1 


-i 1 



&2 



o\- L 

1 1 



i 1 



-i 1 
i 1 



-i 1 



V 



i 1 



1 1 .. 



-i 1 
..0 



By interchanging the rows, we obtain 



li 



detRo 



-l) M det 



/ -i 1 
i 1 



V 



-1 



-i 1 
i 1) 



(-1) 



M 



1 



22M+1 



i 2M+1 . 



(3.18) 



It follows from (3.14), (3.16), (3.17) and (3.18) that 

det R p = det(i? 3 ) det(J> F P ) det R? det i?f = -1. 

The calculation of det R A is similar. The linear transformation R A is a composition of 
transformations, 

R A . C 2{2M+^ A ^IX C 2(2M+1) C 2(2M+1) C 2(2M+1) X J^?c2(2M+l) . 

Here Fa is the finite inverse Fourier transform in the anti-periodic representation, R A is 
the transformation from the Fourier series representation to {e A k , e A _ k } k ^ = _ M , while R A 
is the transformation from {a A * ,a A k }^ = _ M to {q(zA{k)),p(zA(k))} k L_ M . The matrix of 
T A ■■ C 2M+1 -> C 2M+1 is given by 



Fa ■-- 



V2M + 1 



M 
A,-M 



y -M 
'AM 



\ 



\ Z A,~M 



,M / 



2ni(k+^)l 

where z l Ak := z l A {k) = e 2M+1 f or = —M,...,M. We recognize this matrix as the 
Vandermonde matrix. It follows from [15] that 



det Ta 



vmJ - 1 2M+l(z A _ M ....^ AtM 
1 



—M 



...Z, 



M \2M+1 



1 Z 



1 Jl 
A-M Z A-M 



1 Z A,M 



-A,M 



y 2M 
"A-M 



y 2M 
-A,M 



("1 



V 2M + 1 -M<j<k<M 



Y[ ( W k+ ^-w j+ h, 



where w = e 2M + 1 . The right hand side of the equation in (3.19) can be written 



(-i) 



M 



V2M+1 



1 ZM+i ( — 1 



,M 



n (^-^ 

-M<j<fe<M 



(3.19) 



(3.20) 



where we used that 

n = (-D ( ™ +1) = (-D 

-M<j<k<M 

Since J-p clearly also is a Vandermonde matrix, we have 

l 



AY 



det Fp 



V2M+l 2M+1 



n 



(3.21) 



-M<j<k<M 
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It follows from (3.15), (3.19), (3.20) and (3.21) that 

detJU = detJ^p = {-i) M . 

Hence, 

det{F A @F A ) = {-l) M . 

Define for k = -M, ...,M, 



(3.22) 



Ri 



a(z A (k)) 



H > k ' ^2 V iz A {k)a{z A {k))- 1 
Then the matrix of R A is given by 



a(z A (k)) 
-i2A(/c)a(z A (/c))~ 1 



R? 



M 



-M 



where 

det R A 
Furthermore, we have 



M 



H) 



■\2M+1 



n z A(k) 



-M 



-1 



-1 



det R% 



A - det B% 



(-1) 



H) 



•2M+1 



2M+1 



(3.23) 



(3.24) 



22M+1 

Combining (3.14), (3.22), (3.23) and (3.24) we obtain 

det R A = det(i? 3 ) det{F A J^a) det R A det i# = 1, 

and the theorem is proved. □ 

Let A and Op denote the unit vacuum vectors in Alt(W+) and Alt(W+) respectively. 
We want to choose a representation of the transfer matrix V A in the Fock representation 
Alt(Wjj5) such that the vacuum vector 0,4 is an eigenvector for V A associated with its largest 
eigenvalue. Similarly, we want to choose a representation of the transfer matrix V p in 
the Fock representation Alt(WT) such that the vacuum vector Op is an eigenvector for V p 
corresponding to its largest eigenvalue. By Theorem 3.1, we can write V as the map, 



V J 



V P : Alt eV en(^+)eAlt 

even 



where 



V A ~ V\ 



Alto 



(Wf) -> Alt eV en(W. 
and V P ~ V 



Alt even (I<), 



Alteven(VF^)" 



Let denote the restriction of T A to Wjj5, and let T p denote the restriction of T p to 
W p . Define as in [13] the linear transformations, T(T A ) and T(T P ), acting on Alt(W+) and 
Alt(IV^) respectively as 

(TfgTf)©...® (T A ® • 



r(Tf ) = i 



T A ■ 



^) 



(3.25) 



2M+1 



and 



r(T p ) = l e t p e (rf ® rf ) 



(^ 



(3.26) 



2M+1 



It can be checked that 



r(r(7£)) = T' 



and T(r(T^)) = T . 
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It follows that there exists two real numbers, A^ and A^, such that the representations of V A 
and V p in the Fock representations Alt(W+) and Alt(Wjf) respectively are given by 

V A = \£r(T?)\ Mteven{wt) and V p = X p T(T p )\ Mtcvcn{wf) . (3.27) 

Since the spectra of the induced rotations T A and T p do not contain 1, both T A and T p 
are strict contractions. Then it follows from (3.25) that the largest eigenvalue of T{T A ) is 1 
with the unique eigenvector given by the vacuum vector Oa in Alt(W+). Hence, A^ is the 
largest eigenvalue of V A with corresponding eigenvector Oa- A similar argument shows that 
A^ is the largest eigenvalue of V p with corresponding eigenvector Op. We now determine the 
eigenvalues for the transfer matrices V A and V p '. We start by proving the following lemma. 

Lemma 3.2. We have 

det(r(Tf)) = (detTf) 22M and det{T{T A )) = {detT A f M for M > 0, (3.28) 

where T(T A ) andT{T P ) are defined in (3.25) and (3.26). 

Proof. We prove the first equation. Let us consider the tensor product (T_T)® fc which acts on 
A\t k {W^). We have that dim(Alt fc (W_f )) = { 2M k +1 ) for k = 1,...,2M + 1. The eigenvalues 
of (T P )® k is a product of e ~^ Zp( - h ^~^ Zp ^-~^ Zp( - lk ^ with a certain multiplicity, where 
— M < h < l 2 < ■ ■ ■ < h < M. Fix a choice e~ J ( Zp ( n ^ in this product. We want to determine 
the multiplicity of e _7 ^ p ^ n ^. That amounts to figure out the number of ways we can combine 
e -i{zp{n)) w jth the remaining factors e~^ Zp ®fi for I = —M,...,M and I ^ n in the product 
above. This is the same as the number of ways we can pick out (k — 1) elements of e~ J ^ Zp ^ 
from 2M elements, which is (f^) • We obtain 

/ 2M \ 

det(T p ® ■ ■ • ® T p ) = ^ e -^ P (-M)) e -j( ZP (-M+i)) . . . e - 7 (*p(JiO)^ ^ = (det(T p )f™). 
k 

It follows that 



r-> n v^2M (2M\ n „2M 

det(T(T p )) = [det(T p )f^( 3 ) = (det(Tf)) 22M , 
and hence the lemma is proved. □ 

Recall that and £p denote the anti-periodic and periodic spectrum respectively. 
Theorem 3.3. The eigenvalues of the transfer matrices V A and V p are given by 

exp [i(± 7 (z A (-M)) ± .... ± l(z A (M)))] (3.29) 

and 

exp [|(± 7 (z P (-M)) ± .... ± j(z P (M)))] (3.30) 

respectively, where 7 is the positive root of 

cosh 7(2) = c\c\ — s* 2 si z+ ^ for z£Y>a or z £ Sp. 

There is an even number of minus signs in both spectra for T < Tq ■ The largest eigenvalues 
of the transfer matrices V A and V p are given by 

A TT p TT iM 

A£ = ] I e 2 and Aq = | [ e 2 

respective/?/. 
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Proof. Recall that the transfer matrices V A and V p given in (2.6) and (2.7) have factors of 
the form Y\ e pq - We have 

dete Ew = e ^ tr{pq) . 

Using the Clifford relations and properties of trace, we have 

tr{pq) = tr(qp) = - tr(pq) 

which implies that tr(pq) = 0. Hence det e^ pq = 1 and it follows that det V A = det V p = 1. 
From (3.27) we have the identity 

detU A = det[A^r(T^)]. (3.31) 
Since A\t(W A ) has dimension 2 2M+1 , and T{T A ) acts on Alt(W^), it follows from (3.31) that 

. 9 2M+1 

1 = (Aq ) det(r(2^)). (3.32) 

From Lemma 3.2, we have 

detr(r A ) = (detT^) 22M , 
and combining this with (3.32), we obtain 

X A _ 1 

Since the eigenvalues of Tjj 4 are given by the set {e _7 ^)} z6 E A , we obtain 

Ao = n e 2 ■ 

It follows from (3.27) that an eigenvector for V A is on the form, e A j 1 A e+j 2 A ... A e A ■ , 
for — M < ji < j2 < • • • < jk < M, where k is even. The corresponding eigenvalue is 

exp ( ^ ^ 7(2;) — y^ y j(zA(ji))j ■ Hence (3.29) follows. The calculation of the eigenvalues 
of y p is done in the exact same way. □ 

4. Spin Operator 

In this section, we compute the induced rotation associated with the spin operator. We 
determine the matrix of this rotation in the orthonormal basis of eigenvectors for T(V). The 
spin operator Oj at site (j, 0) acts on the Clifford generators as 

VjPkvJ 1 = ~ sgn(/c - j - l)p k , 
VjqkVj 1 = - sgn(fc - j)q k , 

where 



, , f +1 if x > 
sgn(x) = | _ x . fx<0 _ 



Since the spin operator <jj anti-commutes with U, it maps the —1 eigenspace for U into the 
+1 eigenspace for U. Thus, we have 



„f A 

of 
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acting on (U = — 1) © (U = +1). Since (U = — 1) © (U = +1) is also unitarily equivalent to 
Alt(W£) and Alt(W£), we can consider aj as the map, 

Uj : Alt(W^) -»• Alt(W^) (4.1) 

or 

o-j : Alt(W^) -> Alt(Wf ). (4.2) 

It is convenient to let the induced rotation for the spin operator be the identity at the 
endpoint j = M. In order to do this, we translate the q k basis elements by 1 lattice unit. 
This will change the operator at site j to be multiplication by — sgn(/c — j — 1) for both q k and 
p k . In the Fourier representation, the corresponding action on the induced rotation for the 
spin operator is multiplication by za on the left and by Zp 1 on the right. The spin operator 
is then completely reduced to a "difference" of periodic and anti-periodic translations. We 
therefore formulate the matrix representation of the induced rotation associated with o~m- 
Let s := T((Jm) denote multiplication by — sgn(/c — M — 1). In the Fourier representation we 
have the following matrix representation for s, 

sf{zA) = y 2 (*az)- m ( } (43) 

J y ' 2M + 1 ^ z - z A 
The induced rotation for the transfer matrix is then given by 

nv)-.= (l ;)too(V J). 

The conjugation of T(V) was not introduced in Section 2, since it complicates the repre- 
sentation of the transfer matrix in the alternating tensor algebra. The eigenvectors of T(V) 
corresponding to the eigenvalues e -7( ^ p ) and e 7<Zp ) are given by 



~ e P (z) ._ J_ ( a(z) 
+ ' k{ ' '~ V2\ Hz)' 1 

and 



8z P {k)(z) (4-4) 



respectively for z £ Sp. The eigenvectors, k (z), corresponding to the eigenvalues, e^ 7 ^), 
are defined in the exact same way for z € XU- We now consider the isotropic splittings, 

VF = W"_f © W p and = © W A , 



where {e^ fc }^f = _ M is the choice of basis for W±, and {e± fc }^f = _ M is the choice of basis for 



W±. The following theorem follows. 



Theorem 4.1. Suppose that 



A B 
C D 



is the matrix of the identity map on W going from the © W p splitting to the © W J _ 
splitting. Then the matrix elements of A,B,C and D are given by, 

1 (z P (k)z A {l))-* 
(2M + 1) z P (k)-z A (l) 



A lk = D lk = — j- - (g^d(0L^ [ o ( ZA (f)) o ( Zp ( fc ))-i + a (z A {l))- l a(zp(k))}, (4.6) 
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Blk = (2M + 1) ^^(OK 1 ^)) " a( ZA (l)Hz P (km (4-7) 

Q* = -B ik (4.8) 

/or fc,/ = —M, ....,M with 



/ A 1 {z)A 2 {z) . 

a(z) = fiF)^P) A = ^7^, j = 1.2 

/or z G Sa o?" 2 G Sp. 

Here we have for example D = Q^l^p which maps W p into W A . Recall that the rep- 
resentations {qkiPk} are generators of an irreducible ^representation of the Clifford algebra 
Cliff (W) on C 2 ( 2M+1 ) [13]. Let F denote the Fock representation associated with these Clif- 
ford relations. Since aMxajJ = sx for x G W, where gm '■ 

C 2(2M+1) ^ C 2(2M+1) ; it foUows 

that 

a M F{x) = F(s(x))a M . (4.9) 

The Fock representation defined in (3.6) is an irreducible ^representation of the Clifford alge- 
bra. All irreducible * representations of the Clifford algebra over W are unitarily equivalent 
[1]. Thus, it follows that there exists a unitary map, 

Up : C 2 ( 2M+1 ) -> Alt(Wf ) 

such that 

and a unitary map, 
such that 



F p (x) = UpF^Up 1 (4.10) 
U A : C5 2M+1 ) -> A\t(W A ) 



F A {s{x)) = U A F{s{x))U A \ (4.11) 
It follows from (4.9), (4.10) and (4.11) that 

U A a M Up 1 F p {x) = F A {s{x))U A a M U p l . 

The map, 

a ■= U A a M U p l : Alt(Wf ) -> Alt(W^), 
is an intertwining map which is unique up to a constant by Schur's lemma. 
Thus, we have the following commutative diagram, 

A\t(W p ) — ^->- Alt(w£) 

F A os 

Alt(Wf)— ^Alt(W-f), 

where a satisfies the intertwining relation, 

o-F p (x) = F A (s(x))a. (4.12) 

We will use this relation in Section 6 when we find an expression for the spin matrix elements. 
Now we can regard a AP as the restriction of aj to Alt cveD (W p ). In a similar fashion, we can 
restrict aj to Alt cvcn (W A ) and recover o PA . This restriction is central since we will consider 
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the spin operator in terms of an eigenvector basis for the transfer matrix V which is simple 
only in Alt eve n(^) and A\t even {W p ). We are interested in understanding the N — > oo 
limit of the two point correlation for the spin operator. The eigenvector associated with 
the largest eigenvalue for the transfer matrix will dominate in this limit under the action of 
the spin operator a. In order to compute this limit, we will need to know the spin matrix 
elements for a. As we will see in Section 6, the spin matrix elements can be written in terms 
of the inverse of D in (4.6). 

5. Two-point Correlation Function 
For n > m, we have [8] 

_ tr(a i V n - m a j V ( - 2N+ V- n+m ) 

\ a rni<7nj)A - tr (y2JV+l) " 

Recall that 

V = V A © V p , 

where 

V A = Xq r(T+)| Altcvon(H M), V P = \ P T(T+ )| A lt ev en(W^)> 

and where T(T A ) and T(T P ) are defined in (3.25) and (3.26). Using the equation above, the 
identity T(T a )0a = Oa, and the fact that only the terms which involve the vacuum vector 
Oa survive in the semi-infinite volume limit, N — > oo, we obtain 

lim {(Jmicrnj)A = (0A,o"i(pr) n "VjO^) 

N— >oo 

= J(n-m)J2fL_ M [ 7 (*P(J))-7(*A(J))] x 
M 

K=l -Af<fci<- 
-<k 2K <M 

x e'^^^\el M A ••• A e^ 2 ^,,0 A ). 

On the torus, where both M and N are fixed, the two-point correlation function can be 
written in a form involving the spin matrix elements [e A ^ A - - - A e A t , ae p ki A • • • A e p kj( ) 
for L + K = even. 



6. Spin Matrix Elements on the Finite Periodic Lattice 



In this section, we will provide a formula for the matrix elements of the spin operator in 
a basis of eigenvectors for the transfer matrix V. Recall that after "conjugation" the spin 
operator is a map, 

a : Alt(W^) -> Alt(W^), 

and the induced rotation, s := T(a), associated with the spin operator is the identity on W. 
We defined the matrix of this map to be 



A 
C 



where for example D = Q A \^p. 
bilinear form (-,-), and D~ T : W p 



D 
D 



Wi 



wi 



Let D T denote the transpose of D with respect to the 
->■ W A is the inverse of D T . Here D is invertible when the 
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one-point function {0 A , aOp) is nonzero. Let V denote the collection of subsets of 
{— M, —M + 1, M}. For an element J in V, we write 

J = {Ji, J 2 , Jfc} with Ji < J 2 < ■ ■ ■ < Jk- 

We write # J = A; for the number of elements in J. If R is a 2M x 2M matrix, we let Ri t j 
denote the (#/ + # J) x (#/ + # J) sub matrix of R made from the rows and columns of R 
indexed by I and J respectively. For I,JeV define 



and 



The sets {e^ 7 } and {e^ j} are orthonormal bases for Alt(W±) and Alt(W±) respectively. 
We have the following theorem. 

Theorem 6.1. Let I,JeV and suppose (0 A ,a0p) / 0. Then the spin matrix elements are 
given by 

(eij,ae^j) = (0 A ,aO P )Pf(Rj,j), 
where Pf (Ri : j) is the Pfaffian of the skew symmetric matrix, 



Ri,j 



-D~/ xI D^C JxJ 
The sum is over all such I and J with #L + # J even. 

Proof. Since a satisfies the intertwining relation (4.12), the proof follows from Theorem B.2. 

□ 

Corollary 6.2. For -M < k < j < M and (0 Al a0 P ) / 0, we have 
(0A,ge^ fe e+j) p ! P ! 



(0 A ,o-0 P ) 
(ei k ,ae p 



(e^ k ,BD' 1 e^ J ) = (BD- 1 ) kj , 



+ ,ki u _ / D -r~P \-n-r (61) 

7. Bugrij-Lisovyy Conjecture for the Spin Matrix Elements 
For simplicity, let £p and Yl A now denote the sets | 2M+1 ' ^ e i~ M}^j and 

1 2M+1 ' k e {— M, ...,M} > respectively. We consider the isotropic case and define the 

interaction constant to be /C := )C\ = /C 2 . The function 7(6*) is defined as the positive root of 
the equation 

cosh( 7 (0)) = sinh(2/C) + sinh(2/C)- 1 - cos(0). 
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A. I. Bugrij and O. Lisovyy [2] proposed the following formula for the spin matrix elements 
on the finite periodic lattice for the isotropic case in the orthonormal basis of eigenvectors 
for the transfer matrix, 

( g +,/i A ~ e ih A - A g +,i m > A 1 %M A - A e+,kj 

1 ( . 1 / \ 

m ± v {e h) m' p ~^\ e 3 ) 

= v^n , e n 



1=1 



(2M + 1) sinh 7 (6^) J=i J(2M + 1) sinh 7 (0£.) 



sm o— »- rr sm „ J 



n 0111 i n 

l<i<i'<m S mh ^5 ^- l<j<j'<m' sin h ^_ V_ 

. U 7(^)+7(C) 

smh 1 



II ra^> ( 7 - x ) 



l<i<m,l<j<m' gin - „ fcj 



where 



e = |l-(sinh(2/C))- 4 |i : 



and where the cylindrical parameters £t and v{e ld ) are defined as 
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n, eSA smh(2M) 



u(e") = log 



Here m + m' is even. It can be checked that e ± ^ e •* can be written in the forms, 



V+{e ld ) := e~^~ = 



and 



e-^(A^(e- e ) ~ A(-l)) Uo ! >o^ A fjVv?) - A(e^)) 

\ 7(0) 7(6') ^ j 

\| e-— (A-i(e")-A(l)) n e '>o G E p e"~(A- 1 (^)-A(e^')) 



«(e«) 

y_(e t0 ) :=e — r" 



ejW^) - A-Hl)) IWs P e^We") - A^e^)) 
\J e— (A(fi») - A-i(-l)) n e '>o e E A e— (A(e^) - A-i(e^)) 

where A(e* 6 ') = e 7 ^ for 6 = 6f_ or 9 = 8^, and where we used the short-hand notation 

j (8) := 7(e ie ). Here the square roots are taken with positive real parts. This particular form 
of v(e td ) will play a central role in Section 10 when we consider the factorization of the ratio 
of summability kernels on the spectral curve. 

8. Pfaffian Formalism and the Bugrij-Lisovyy Conjecture 

The Bugrij-Lisovyy formula in (7.1) and Corollary 6.2 lead to the following conjecture for 
the isotropic Ising model: 
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Conjecture. The matrix elements of D T , BD 1 and D l C are given by, 



V + (z A (i)) V-{z P {j)) sinh 



7(0^+7(0 



^(2M + 1) sinh 7 (^ A ) ^(2M + 1) sinh 7 (#f ) sin 



OC 



BD~}, 

V+(z A (i)) V+Mif)) sin_ 



(2M + 1) sinh 7 (6»f ) J(2M + 1) sinh 7 (6^) sinh 



V-(zpU)) V-Mf)) sin _ 

v " " " = ' ' " 7(gf)+7(ep 



'(2M + 1) sinh 7 (0f ) ^(2M + 1) sinh 7 (^) sinh 

respectively, where oc denotes proportional to. 

The proportionality constants in each of the terms above are the same and therefore 
omitted. Here we have extended the Bugrij-Lisovyy conjecture to include the elements 
(e+ k-> a ^+j) wmcn are n °t "physical elements". By applying Theorem 6.1 and Wicks Theo- 
rem, we can express the spin matrix elements (e^ ^ A ... A lm , ae^_ ki A ... A k ,) in terms 
of Pfaffians of a skew symmetric matrix whose elements are multiples of (e^ ^ A f aOp), 
(OA,o~e+ k . Aejj, J and (e+i . , cre^ k . ) . By using the proposed Bugrij-Lisovyy formula for 
these elements, we show that the Pfafnan of the matrix can be written as a product of Jaco- 
bian elliptic functions in the uniformization parameter of the spectral curve. We find that up 
to a constant, this product is the same as the one given in the Bugrij-Lisovyy formula. This 
supports the Bugrij-Lisovyy conjecture. We refer the reader to Appendix C for a description 
of elliptic functions and the Uniformization Theorem C.l. We start by writing the factor 



sinh( 



7(g)+7(Q')- 



sin 



2 > 

in (8.1) in terms of the uniformization parameters by performing the elliptic substitutions, 

e %e = z{u, a) = k sn(u + ia) sn(u — ia), e %e = z{u , a) = k sn(u + ia) sn(u' — ia), 
where 

K' 

%u = Q U ' = — and G (0, 2K), 

and < 2a < K' is defined by s\ = — isn(2za). Here K and K' are elliptic integrals given in 
(C.8) and (C.9), and 
Theorem C.l implies 



(C.8) and (C.9), and sn(u) is a Jacobian elliptic function of u with modulus k = Then 



y+ei 

e 2 = vfcsn(-u — ia), 
e 2 = v ksn(u + ia). 
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We have 

7+7' 7+7' 7+tg 7 / +»6> / 7+i0 j'+id' 



e" — e* 



Thus the numerator of (8.2) becomes 

ik[sn(u — ia) sn{v! — ia) — sn(V + ia) sn(u + ia)]. 

Using the identity 

sn 2 (u) — sn 2 (ia) 

sn(u - ia) snfit + ia) = -5 — — ^— — 

1 — k sn z [u) sn / (^a) 

which follows from the addition formulas (C.4) and (C.5), the denominator becomes 
( , sn 2 (u) — sn 2 (ia) \ / sn 2 (V) — sn 2 (ia) 



1 — k 2 sn 2 (u) sn 2 (ia) / \ 1 — k 2 sn 2 (w') sn 2 (ia) 
After some simplifications, the expression in (8.2) becomes 

— 2i[cn(ia) dn(ia) sn(ia)] [sn(«) cn(u') dn(u') + sn(u') cn(it) dn(it)] 
(sn 2 (u) - sn 2 (n'))(l - fc 2 sn 4 (za)) 

= si , 1 ,y (8.3) 
sn(u — u ) 

In the isotropic case we have s^ 1 = y/k. If we substitute v := — u + and v' := —u' — 
in (8.3), and apply identity (C.10), we obtain 

si — — = — \/ksn(v — v). 

sn(u — w) 

Now we will make the following elliptic substitutions, 

i0 A 

e l i = fcsn(uj. — ia) sn(uj . + ia) for 1 < j < m, 

id p 

e k i = ksn(u kj — ia) sn(u kj + ia) for 1 < j < m , 
followed by the translations, 

iK' 

vij : = -u tj + — for 1 < j < m, 

iK' r . 
v kj : = -u kj — for 1 < j < m . 

Then using Theorem 6.1 and Wicks Theorem, we obtain for m + m! = even, 
<< h A e£, a A . . . A e$ u , ae% M A e+ k2 A ... A e^, ) 
= (0^,<70 P )Pf(r), 

where Pf (r) is the Pfaffian of the skew symmetric matrix r with matrix elements above the 
diagonal given by: For 1 < i < i' < m, 

(ei tli e^ <tOp) 

r kh> = 



(2M + 1) sinh 7 (0£) J(2M + 1) sinh 7 (0/J ) 



V+{z A (li)) V + (z A {h)) , n- f „ 

oc - {-Vksn(v k - vi ,)), 
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for 1 < i < m and 1 < j < rri , 
r k , m+kj ( 0a;Ct0p ) 

oc - mm . mm ( -^ mK - ^» 

yJ(2M + 1) sinh 7 (0jj) yJ(2M + 1) sinh 7 (0£ ) 
and for 1 < j < j' < rri , 

r m+ k 3 , m+kjl - (0A)(t0p) 

a mm p v JsML=(-^ (vtl - % ,)). 

^/PM + l)si„h 7 (»£) + l)sinh 7 (Sp 

Define to be the (m + rri) x (m + rri) diagonal matrix 



F- 6 ° 
E ~ 



with matrix elements 



ef, = ; V+{ZAih)) for !<*<m, 

l * 1 v /(2M + l)sinh 7 (^) 

e^ +fc . m+fe . = V V J77 for 1< j < rri. 

m+k ^ +k] ^(2M + l)sinh 7 (^) 

Define s to be the (m + m') x (m + m') skew symmetric matrix with matrix elements above 
the diagonal given by 

S U;/ = -Vfcsn^ - uj.,) for 1 < i < i' < m, 
sji.m+fcj = -Vfcsn^ - u fej .) for 1 < i < m, 1 < j < rri, 
Sm+k^m+kj, = -Vksn(v kj - v kjl ) for l<j<f< rri . 

We will need the following lemma found on page 87 of [13]. 

Lemma 8.1. [13] Let r be the 2n x 2n skew symmetric matrix with i,j matrix element 
Tij = —\/ r ksn(ui — Uj) for i,j = 1,2, ...,2n. Then 



2n 



Pf (r) = JJ r. 
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Using Lemma 8.1, we have 

Pf(r) = Pi(EsE T ) 
= (det-E)Pf(s) 

_-rr V^Mji)) tt ^(fHM) 

'= i ^(2M + l)sinh 7 (^) }=i ^(2M + l)sinh 7 (^.) 

x JJ - Vk sn(v h - vi., ) Yl sn (% ~ v k 3 , ) x 

l<i<i'<m l<j<j'<m' 

x HI -Vksn(v h - v kj ), 

1< i<m,l<j<m' 

which up to a constant is the same product formula as given in (7.1). 

9. Numerical Calculations 

We have numerically compared the calculation of our BD- 1 , D~ T , and D^C matrix 
elements with the corresponding terms in the Bugrij-Lisovyy formula (7.1). What we find is, 
that as we let the temperature approach Tc from below, our results are close to theirs, but 
there is a discrepancy. We are uncertain of the reason for this. We also find that in the scaling 
regime, i.e. as we let M get bigger while we keep the temperature close to Tq, the precision 
improves. We compared the formulas in MATLAB by creating our (2M+1) x (2M + 1) matrix 
D as given in (4.6). Then we multiplied the transpose of this matrix with the matrix with 
elements (e^ k , ae^ •) for —M < k,j < M as given in the Bugrij-Lisovyy formula (7.1). If this 
had returned the identity matrix, it would have shown numerically that the Bugrij-Lisovyy 
formula provides an inverse for our matrix D T . 

10. Holomorphic Factorization of the Ratio of Summability Kernels on the 

Spectral Curve 

In this section, we exhibit the "new" elements V± in the Bugrij-Lisovyy formula as part 
of a holomorphic factorization of the periodic and anti-periodic summability kernels on the 
spectral curve associated with the induced rotation for the transfer matrix. 

Introduce the cycles 

N + = {m|0 < < 2K,Qu = 0}, 

A/1 = {u\0 < < 2K,%u = ±K'}, 

where K and K' are elliptic integrals given in (C.8) and (C.9). These cycles are the boundaries 
for the holomorphic factorization of the summability kernels. Introduce the notation, A p := 
\{e tp ) = e 1 ^. We will prove the following theorem. 

Theorem 10.1. Near the curve, 

AA_ = [u : < Ru < 2K, <3u = ±K'}, 

we have the identity 

V+(u) V-{u) 
z(u) 2M + 1 + 1 z(u) 2M + 1 - 1 

while near the curve, 

Af + = {u:0<Mu< 2K, = 0} , 
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(10.1) 



we have the identity 

V+{u) V-(u) 



z ( u )2M+l + l Z (n)2M+l_l' 

where 



(10.2) 



V+(u) 



f (A(n)-A^)n p >o e s A (A(n)-A p ) 
(A(«) - A ) n p >oeSp( A («) - a p) 



VlAW-A^^n^oeE^AW-Ap 1 )- 

27ie square roots are here chosen to have positive real parts. 

Before we give the proof we explain our principal concern with Theorem 10.1. Recall from 
Appendix A that the two cycles A4+ and A4_ on the spectral curve Ai, where A < 1 and 
A > 1 respectively, are given by 

M ± = {{z,\) = {e w ,e^)}. 

Recall from Appendix C that in the uniformization parameter u, the cycles M± are located 
at the following positions in the periodic parallelogram, 

M± = {u : < < 2K, 9u = ±4f }. 

We observe that in a neighborhood of the cycle M.+, the function V+(u) is a multiple of 
V+(u) in (7.2) from the Bugrij-Lisovyy formula while V-(u) is a multiple of V-(u) in (7.3) in 
a neighborhood of the cycle M-. The factorization problem solved by V± makes it possible 
to transform sums over the periodic points to the anti-periodic points. A related factorization 
in the scaling limit was employed in [10] to compute a relevant Green function for the Dirac 
operator on the cylinder. Subsequently, in [14] the technique in [10] was used to compute 
the matrix elements of the spin operator for the periodic Ising model in the scaling limit. It 
is expected that this technique can be extended to the finite periodic lattice to provide an 
alternative proof of the Bugrij-Lisovyy conjecture. If one uses formulas for the spin matrix 
elements on a finite periodic lattice, one can control the convergence of the scaling limit in 
a much simpler way. This convergence result is given in [7] assuming the Bugrij-Lisovyy 
conjecture [2]. 

We now return to the proof of Theorem 10.1. 

Proof. Recall from Appendix C that < 2a < K' is defined by S\ = — isn(2ia). A simple 
calculation shows that a substitution of a with a + ^- in the elliptic parametrization of the 
Boltzmann weights, interchanges z and A in the spectral curve (C.l), and sends si to — S2, 
S2 to —si, and c\C2 to — c\C2- We therefore introduce the relation 



. %K' 
la := — — h xa. 



Then 



and 



A' := X(u + , of) = z(u, a) = k sn(n + ia) sn(u — ia) (10.3) 



z' :=z{u + tJ ^,a') = \(u,a) = — / '- (10.4) 
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for k = . Interchanging A and z through the substitution, a4c+ M|- , the spectral curve 



given in (C.l) becomes 



A' + A'- 1 , z' + z'-i 
si \~s 2 =cic 2 . (10.5) 



Using the spectral curve in (10.5), we can write 

| [ (C1C 2 - Sl COSp - 52^1—) = J [ ^l( V+ 2 / ^ - s l cos P 



= (t) 2M+1 II " eiP )( 1 " A' _1 e" ip )) 
peSp 

= ( a ) 2Af+l (A ,2Af+l _ 1)2 y-(2M + l) (1Q 6) 



and similarly, 



J] (ciC2 - -i cosp - s 2 ^) = (f) 2M+1 (A' 2M+1 + 1) 2 A'- (2M+1) . (10.7) 

Then from (10.3), (10.6) and (10.7) we obtain 

' z {u) 2M+1 + \ \ 2 
z (u) 2M + 1 - 1 J 

y2M+l , ^2 



A /2M+1 _ 1 

n pg s A ( c i c 2 -sicosp- S2 z,+ f 

„/-l. 



(10.8) 



n p6 E P ( C l C 2 - glCOSp- S 2 Z ' +Z 2 ' j 

Now we rewrite the factor in (10.8) using (10.4) and the spectral curve in (C.l). It becomes 

rwAp + A-j-^ + o) 



ru Sp (A P + A p - i -(^ + ^- 1 )) 
riWAp-^u-Ap-v- 1 ) 



" ripeEpCAp-^a-ApV- 1 ) 

n pe s A (AH-A p )(A(^)-A- 1 ) 
"n peS p(AH-A P )(A(n)-A p 1 )' 

Since A(p) = A(— p) for p / 0, 7r, the right hand side of the equation above can be written 

(A(u) - A 7r )(A(n) - A^ 1 ) E[p>o e E A (AH - A p ) 2 (A(u) - Ap 1 ) 2 



Now define 



(A(n) - A )(A(«) - Aq 1 ) n p >o e Ep(AH - A P ) 2 (A(n) - Ap 1 )^ 
V+(u) :-- 



'(\(u)-K)Up>oev A M u ) ~ A P ) 



(AH - A ) rip>06Ep(AH - A p ) ' 



, (AH-A 1 )n p>0g Sp(AH-A p - 1 ) 
(AH-A^^n^oeE^AW-Ap 1 ) 
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such that we have 

) 2 = ( YzM V (10 9) 

In order to figure out the correct signs near the curves J\f± when we take the square root in 
(10.9), it is enough to check the identities in (10.1) and (10.2) at one point. We will check 
(10.1) by letting u = K ± iK' which are points on the curve JV_. Using the identities (C.10), 
the addition formula (C.4), and Theorem C.l, we obtain 

z(K ± iK') = k sn(K ± iK' + ia) sn(K ± iK' - ia) 

, ns(-fT) cn(ia) dn(ia) — A: -1 ds(K) cs(K) sn(ia) 



1 — ns 2 (K) sn 2 (ia) 
ns(if) cn(za) dn(za) + k~ l ds(K) cs(K) sn(ia) 
1 — ns 2 (K) sn 2 (ia) 

Here we also used the fact that sn(it) is an odd function of u while cn(u) and dn(it) are even 
functions of u, [17, p. 493]. Using the fact that sn(K) = 1, cn(K) = [17, p. 499], the 
identities in (C.2) and the fact that «2 = — fc~ 1 ns 2 (ia) [13, p. 338], we obtain 

z(K ±iK) = - j— 

1 — sn^(za) 

_ fc-^l - A; 2 sn 2 (m)) 

1 — sn 2 (za) 

_ fc _1 (H- ka^ 1 ) 
1 + /c _1 «2 1 

which can be checked to be greater than one for T < Tq. Thus, z ^±iK^) iM+1 -i * s P os itive. 
A similar calculation gives, X(K ± iK') = 1, and hence V+(u) and V-(u) are both positive 
for u = K ± iK' . Thus the identity in (10.1) follows. Now we prove (10.2) by letting u = K 
which is a point on the curve A/+. Then it can be checked that 

kcn 2 (ia) , .„,._] 

< K ) = A 2 f - ; = ± zK ) < 

dn (ia) 

and hence z( - E -- ) ali + i_ 1 is negative. We also have \(K) = 1 so that V±(u) are positive for 
u = K. Hence (10.2) follows, and the theorem is proved. 

□ 

The author would like to thank Professor John Palmer for his advice and support of this 
work. 

Appendix A. Grassmann Algebra and Fock Representations of the Clifford 

Algebra 

In this section, we introduce the Fock representations of the Clifford algebra. We follow 
[13] closely and refer the reader to this work for more details. Assume W is a finite even- 
dimensional complex vector space with a distinguished nondegenerate complex bilinear form 
(•,•). A subspace V of W is isotropic if (x,y) = for all x,y £ V. We are interested in a 
decomposing of the space W into two isotropic subspaces W± such that 

W = W + @W^. (A.l) 
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Such a splitting is called an isotropic splitting or a polarization. We parametrize each isotropic 
splitting by an operator Q defined by 

Qx={ X !! X G (A-2) 
^ \ -x if x G W-. 1 ; 

Define Q± := \{I ± Q). Since Q 2 = /, we observe that Q± are the projections onto the ±1 
eigenspaces for Q. We notice that Q + + Q~ = I and Q+Q- = 0. If we define W± := Q±W, 
then the space Fy is the direct sum given in (A.l). An operator Q that is skew symmetric 
and satisfies Q 2 = I is called a polarization. Let Sk denote the group of permutations on the 
set {1, 2, k}. The linear operator defined by 

W m 3 w ' ^ alt(if) := — V] sgn(cr)tt; (Jl (g> • • • <2> 

is a projection from W® k onto Alt fc (VF), where Alt fc (VF) is the space of alternating k tensors 
over W. The wedge product of v G Alt k (W) and w G Alt l (W) is here defined by 

v A w := alt(v ® w) G Alt k+l (W) 

\Jk\yl\ 

and it follows that for Vi £ W and i = 1, k, 

vi A w 2 A • • ■ A v k = —r= ^2 s § n (°')^i <S> v U2 (g> ■ ■ ■ <g> 

(see [13] and [16]). The Clifford algebra Cliff (W) of W is the associative algebra with multi- 
plicative unit e, generated by elements x £ W that satisfy the Clifford relations, 

xy + yx = (x, y)e for x,y £ W. (A. 3) 

For each polarization Q of the isotropic splitting of the space W, there is a Fock representation 
Fq of the Clifford algebra Cliff (W). This representation acts on the alternating tensor 
algebra, 



Alt(W+) :=0Alt fc (^ + ), 



k=0 

where Alt°(W+) = C and n = dim(W+). The Fock representation is defined as 

W 3 x H> Fq(x) := c(x + ) + a(x-) 

for x± = Q±x. Here W- is identified with the dual Wl via the nondegenerate complex 
bilinear form W + 3 x + t— > (x + ,x-) for x_ G The creation operator c(x+) associated 

with x + £ W+ acts on Alt fc (W+) in the following way, 

k\t k {W + ) 3»4 c(x+)w = i + Ad£ Alt fc+1 (M/+). (A.4) 

The annihilation operator a(x_) associated with x_ G VF_ is defined as a(x_) := c T (x_), 
where c T (x_) is the transpose of c(x_) with respect to the complex bilinear form (•,•)• It is 
given by 

k 

a{x-)v = ^(— l) J ' -1 (a;_, Vj)v\ A ■ ■ ■ Avj A - ■ ■ Avk (A.5) 
j'=i 
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for v = v± A • ■ ■ A Vj A ■ ■ • A G Alt fc (W + ), and where Vj signifies that the factor Vj is 
omitted from v. It can be checked that the creation and annihilation operators satisfy the 
anti-commutation relations, 

c ( x +) c (y+) + c (y+) c ( x +) = o, 

a(x_)a(y_) + a(y_)a(x_) = 0, (A. 6) 

a(x-)c(y + ) + c{y + )a(x_) = (x-,y+)I 

for x±,y± G W± . Since W± are isotropic subspaces and by using the anti-commutation 
relations given in (A. 6), it is not hard to see that Fq satisfies the generator relations for the 
Clifford algebra, 

F Q (x)F Q {y) + F Q {y)F Q {x) = (x,y)I for x,y £ W. (A.7) 

When the polarization is understood, we will drop the subscript Q and write F := Fq. We 
write 

o := leoe- • -eo 

for the vacuum vector in Alt(W+). The vacuum vector is defined to be the unique vector 
that is annihilated by all the elements in W- in the Fq representation of the Clifford algebra. 

Define the Clifford group Q as the group of invertible elements g in the Clifford algebra 
Cliff(W) that satisfy 

gvg- 1 = T{g)v for v G W C Cliff(W) 

for some linear map T{g) on W. It follows from this equation that T is complex orthogonal, 
i.e. 

(T(g)v,T(g)w) = (v,w) for v, w G W. 

Appendix B. Berezin Integral Representation for the Matrix Elements 

In this section we introduce a representation of the creation and annihilation operators 
which is an analog of the holomorphic representations as given in Faddeev and Slavnov [4] . 
We will use this representation to write the matrix elements for the Fock representation of an 
element g in the Clifford group as Pfaffians of a skew symmetric matrix whose entries depend 
on the matrix elements of the induced rotation associated with g. 

Assume W is a finite even-dimensional complex vector space with a Hermitian inner prod- 
uct (-, •) and a distinguished nondegenerate complex bilinear form (-, ■) defined by 

(«, v) = (u, v) for u, v G W, 

where u H> u is a conjugation. We consider a Hermitian polarization, 

W = W+®W^, 

where W± are isotropic subspaces of W as defined in Appendix A. Let {e^"} and {e^} denote 
orthonormal bases for W+ and W- with respect to the Hermitian symmetric inner product, 
and suppose that {e^~} and {e^} are dual with respect to the complex bilinear form, i.e. 
( e fe > e i~) = ^ kl ' Suppose that W has dimension 2M and define 

ef := e J A • ■ ■ A e J for 1 < h < ■ ■ ■ < I k < M. (B.l) 

The set {e^} is then an orthonormal basis for Alt(W±), where we define = 1. Let V 
denote the collection of subsets of {1, ....,M}. For an element J in V, we write 

J = { Ji, J 2 , Jk] with Ji < J 2 < ■ ■ ■ < Jk- 
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We write # J = k for the number of elements in J. If R is a 2M x 2M matrix, we let Ri t j 
denote the (#/ + # J) x (#/ + # J) submatrix of R made from the rows and columns of R 
indexed by / and J respectively. An element in Alt(W+) is given by 

G(e+) := £ G je + 
lev 

where the map, V 3 I — > G[ G C. Here G(e + ) can be thought of as a "polynomial" in 
the elements ef in the exterior algebra. For 1 < i < M, the creation operator c(-) acts on 
Alt fc (IF+) as 

Alt fc (IF+) c(ef)v = ef Ave Alt k+1 {W + ) 

for ef G W+. For 1 < i < M, the annihilation operator, a(ef) := -^r , is analogous to a 
"derivative" , and is the linear map 

^ : Alt(W+) -+ Alt(W + ) 

defined as follows: If the monomial, A := e+ A e + A • • • A ef , contains exactly one factor ef 
then 

d 

— j-(ef A ef A ■ ■ ■ A ef ) = ±ef A • • • A ef A • • • A ef , 

where ef signifies that the factor ef is omitted from X, and the plus or minus sign is 
determined by number of interchanges the operator has to make from the left before it 

contracts the factor ef. An even number of interchanges gives a positive sign and an odd 
number gives a minus sign. If the monomial does not contain any factor ef or if it contains 
more than one factor of ef, then X = 0. We define Berezin integrals as linear functionals 

in the following way, 

/.♦*♦ = !, /.-*- = i, jV = o, /*- = o, 

where we assume that de~ and de + anti-commute with each other as well as with and e + . 
An element in the Grassmann algebra Alt (IF) is given by 

G(e + , e~ ) = G 00 + Goie^ + Gioe^ + G u e^ A ef + ... + G 1: ... :M ,M,-A e i A .. A e M A A ... A e^, 

where Goo, Goi, Gio and Gi J ... ) Af,M,...,i are complex numbers [4]. The integral of G(e + , e~) 
is then defined as 

M 

G{e + , e~) Y\ defde^ = G; m,m...i- 

k=i 

For G and H in Alt(W+), we have [4, p. 53] 

. M 

(G, H) = / G(e + )H{e + )e-^ e t^\[defdel, (B.2) 
^ k=i 



where 

G(e + ) = ^G /e 7 

and the conjugation is defined as 



G/eJf A • ■ ■ A ef = G I e 1 A ■ • ■ A e M . 
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It is here understood that ^ e t he k is the power series in the exterior algebra, 

M 

E 



k\ 

k=0 x 

We are in particular interested in two polarizations, 

W = W A © W A and W = W p © VF P . 
Here W± and PV^ are isotropic subspaces of W defined by 

w£ = Q A W and = Q^W, 

where 

Q A :=±(I±Q A ), Q P :=\(I±Q P ), 

and where Q A and Q p are polarizations as defined in (A. 2). Let F p and F A denote the 
Fock representations associated with the Clifford algebra Cliff(W) acting on Alt(W+) and 
Alt(W+) respectively. We consider a map 

g : Alt(W p ) -»• Alt(W^), 

which satisfies the intertwining relation 

gF p {x) = F A {T{g)x)g 

for x € W, and where T := T(g) is the induced rotation associated with g. We write 



T{9) :-- 



A B 
C D 



for the matrix of the induced rotation for g, where T is a map W p © W p — > W A © WL 
Since T(g) is a complex orthogonal matrix, we have the relation 

T(gYT(g) = I, 

where 

D T B T 



This relation implies the following identities, 



D T A + B T C T = I, D T B + B T D = 
C T A + A T C = 0, C T B + A T D = I 



(B.3) 



when D is invertible. Introduce the anti-commuting "dummy" variables af £ W p which 
anti-commute with ef as well. The action of the operator g on G(a + ) G Alt(W+) is given 
by [4, p. 55] 

. M 

( g G)(e + ) = I g{e + ,a-)G{a + )e-^=i a * Aa * Y[da+da^. 
J k=i 

Let {ef } and {af} denote orthonormal bases for W± and W± respectively and for I,JeV 
define ef and af as in (B.l). We first prove the following lemma. 
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Lemma B.l. Let g(e + ,a ) := e^, where 

n = ^Z\(k a im e i A e m) + ( b i m e-t A a m ) + {\ci m aj A a m )] 



Lm 



A B 
C D 



and a = BD 1 , b = D T and c = D l C ' , where 

T = 

is the matrix of a complex orthogonal map, W+ © W_ ->■ © W_. Then 

r M 

gG{e + )= / g{e + ,a~)G(a + )e-^i a k Aa * Y[da+da k 
-* k=i 

defines a linear map which satisfies the intertwining relation,: 

gF p (x)G = F A (Tx)gG 

for G e Alt(Wf ) and x £ W. 
Proof. For af,G€ Alt(W£), we have 

(gF p (al)G)(e + ) = j g(e + , a')a+G(a + )e- ^f=i <*t*°i f[ da+da' 



M 



k=l 



By the signed Leibniz rule, we have 

d 



da. 



e - E2Li = a + e - Ef=i <Aa- 



and it follows that the integral in (B.6) can be written 



f 8 _ M 

{gF p {ajG)){e+) = - J ^— r [g(e + ,a')]G(a + )e~^ a t^k J^da^da 

i k=l 



where 







da. 



^ l i ' 

= (D- T + D~ 1 C)afg(e + ,a~). 



M 



Now we have for v G Alt(W+), 

F A (Taf)v = £ 

It follows that 



k=i 







V. 



(F A (Tat)gG)(e 

M 



Akie~l + Cki , 
L fc=i ae fc 



fc=i 
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We have 

d 



g(e + ,a ) 



= (A + D T B T C + D- 1 C)afg(e + , a~) 

= (D~ T + D- 1 C)a+g(e + ,a-), 
where we in the last equation used (B.3). Thus, we have showed that 

gF p (af)G = F A (Taf)gG 
for af, G G A\t(W p ). We can prove (B.5) in a similar fashion for G W p . □ 

Let Op and Oa denote the unit vacuum vectors in Alt(W+) and Alt(W+) respectively. 
Since the map g satisfies the relation in (B.5), the range of g is invariant under the action of 
the Fock representation F A (x) for all x G W. The Fock representation is irreducible [13], so 
the range of g is either trivial or all of Alt(W+ ). Since gOp is nonzero, the range of g must 
be Alt(W^5). Thus g is invertible, so we have that g is an element of the Clifford group. 

Recall (see [13]) that the Pfaffian of a 2M x 2M skew symmetric matrix R with matrix 
elements Rj^ is defined in the following way: Let {ej} denote the standard basis of C 2n . The 
Pfaffian, Pr(ii), of R is defined by 

^ / 2M \ M 

oWTn ( E R ^ Aek ) = Pf ( i? ) ei A • • • A e2M - ( BJ ) 
• K j,k=i J 

We use Lemma B.l to prove that g can be written as an exponential, whose argument is a 
quadratic form. 



Theorem B.2. Suppose that g satisfies the intertwining relation, 



gF p {x)v = F A (Tx)gv for xeW and v G A\t(W^ 



where 

is the matrix of its induced rotation, T : W p © W p — > W A © W A . Suppose that the one-point 
function (Oa, gOp) is nonzero. Then the kernel g(e + ,a~) of g can be written as 

g(e + ,a') = (0 A ,gOp) £ Pi(Rr t j)ef A ctj, 
i,Jev 

where Pf(i?/ t j) is the Pfaffian of the (#/ + jf J) x (#/ + # J) skew symmetric matrix, 

n D iL)- ^ 

The sum is over all such I and J with + jfJ even. 
Proof. Since g satisfies the intertwining relation, 

gF p (x) = F A (Tx)g for x G W, 
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Lemma B.l implies that there is a nonzero constant A such that 

g(e + ,a') = Xe n , (B.9) 

where 

l,m 

for a = BD^ 1 , b = D~ T and c = D~ 1 C. It can be checked that 

r 1 M 

fc=i 

Writing out the Taylor series expansion in the integrand expression above, we obtain 

(0 A ,~g0 P ) = l. 



Hence, 
and 

It is well-known that 



A = <0 A ,sOp> 

5 (e + ,cO = (0 A ,gO P }e n . (B.10) 



e^= J] Pf^e+AaJ, (B.ll) 
J,Je7> 

where 22j j is given in (B.8). Here + # J must be even in order to contribute to the sum. 
Combining (B.10) and (B.ll), we obtain the result. □ 

Appendix C. The Spectral Curve associated with the Induced Rotation for 

The Transfer Matrix 

In this section, we recall facts about the Jacobian elliptic functions, sn(iz, k), cn(u,k) and 
dn(u,k), where u is the uniformization parameter and k is the modulus. These functions 
play a central role in the holomorphic factorization of the ratio of summability kernels on the 
spectral curve in Section 10. They are also a key element in a product formula for the spin 
matrix elements on the finite periodic lattice as we will discover in Section 8. We follow the 
introduction of the Jacobian elliptic functions as given in [17] and [13], and refer the reader 
to these books for more details. The spectral curve M. associated with the induced rotation 
T Z (V) for the transfer matrix is given by the set of (z, A) such that 

z + z' 1 A + A- 1 
si hs 2 =cic 2 . (C.l) 

The spectral curve is topologically a torus and the two fold covering, (A, z) >->■ z, is ramified 
at z = af,af, where [13, p. 65] 

a.\ = (c* - s*)(c 2 + s 2 ) and a 2 = (c* + s^)(c 2 + s 2 ). 

The roots ol\ and a 2 were introduced in Section 2 in connection with the Boltzmann weights. 
We use the shorthand notation sn(u), cn(u) and dn(n), for sn(u, k), cn(u,k) and dn(u, k) 
since the modulus k is fixed at k = — !— in our calculations. The functions sn(V), cn (it) and 
dn(it) are doubly periodic meromorphic functions of it and they satisfy the equations 

sn 2 (it) + cn 2 (it) = 1, dn 2 (u) + k 2 sn 2 (it) = 1 (C.2) 
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and 

d 



sn(u) = cn(u) dn(u). (C3) 
du 



We use the standard notation [13], 



1 , s cn(u) 

ns[u) := — cs(u) :- 



sn(u) ' sn(u) 
and in general 

nx(u) := — and xy(w) := — ^—l, 

1 ; xn(u) Jy ' yn(«)' 

where x and y are one of either c, d or s. The Jacobian elliptic functions satisfy the following 

addition formulas [17, p. 496-497], 

m r„ i „^ _ sn(") cnfo) dn(v) + sn(u) cn(u) dn(ix) 

SU{U -+- V — ~2 r J--- , l^- 4 J 

1 — k z sn z (u) sn z {v) 
sn 2 (u) - sn 2 (t>) 

sn(u-'u) = — — -— - ; . ' —— — , (C.5) 

sn(u) cn(i>) dn(i>) + sn(v) cn(u) dn(u) 

cn(u) cn(-u) — sn(it) sn(v) dn(u) dn(v ) 
cn(u + = ~2 — 2TI — TT\ ' (^-") 

, x _ dn(u) dnfo) - fc 2 sn(n) sn(t>) cn(n) cnfo) 
dI H u + ^- l-fe 2 sn2(n)sn2(t,) " ^ 

Introduce the elliptic integrals [17, p. 501], 

K= C {l-t 2 )-\{l-k 2 t 2 )'h dt, (C.8) 
J 

i^'= / (1 - £ 2 )-l(l - fc /2 t 2 )-2 (C.9) 
JO 

for which the complementary modulus k' is defined by k 2 + k' 2 = 1. If we translate the 
Jacobian elliptic functions by iK', we have [17, p. 503] 

sn(u ± iK') = fc~ 1 ns(u), 

cn\u±iK') =Tik~ 1 ds(u), (CIO) 
dn(u ± iK') = =Fi cs(u). 

Introduce the two cycles on A4 , 

M± = {(z,\) = (e M ,e^)}, 

for 6 G [— 7r,7r). Here we have introduced the notation 7(0) := 7(e* e ), where 7 is defined as 
the positive root of 

cosh 7(2) = c* 2 c\ — s 2 s i( z+ 2 

Palmer [13, p. 71] showed that an elliptic substitution gives a uniformization of the whole 
complex curve: 

Theorem C.l. [13] The map, 

[0, 2K] x i[—K', K'] 9«4 (z(u, a), A(u, a)) 

a) = k sn(u + ia) sn(u — ia) 
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and 




sn(u — ia) 



sn(u + ia) 



+ S 2 



x + x- 1 

2 



= ClC 2 , 



where k = and < 2a < K' is defined by 

s\ = —isn(2ia). 

In the uniformization parameter u, the cycles Ai± are located at 




References 



[1] R. Brauer and H. Weyl. Spinors in n-dimensions. Amer. J. Math., 57:425-449, 1935. 

[2] A. I. Bugrij and O. O. Lisovyy. Spin matrix elements in 2D Ising model on the finite lattice. Phys. Lett. 
A, 319(3-4):390-394, 2003. 

[3] A. I. Bugrij and O. O. Lisovyy. Correlation Function of the Two-Dimensional Ising Model on a Finite 
Lattice. II. Teoret. Mat. Fiz., 140(1):113-127, 2004. 

[4] L. D. Faddeev and A. A. Slavnov. Gauge Fields, Introduction to Quantum Theory. The Ben- 
jamin/ Cummings Publishing Company, INC, Advanced Book Program, Reading, Massachusetts, 1980. 

[5] G. von Gehlen, N. Iorgov, S. Pakuliak, V. Shadura, and Yu Tykhyy. Form-factors in the Baxter- Bazhanov- 
Stroganov model I: Norms and matrix elements. J. Phys. A: Math. Theor., 40:14117, 2007. 

[6] G. von Gehlen, N. Iorgov, S. Pakuliak, V. Shadura, and Yu Tykhyy. Form-factors in the Baxter- Bazhanov- 
Stroganov model II: Ising model on the finite lattice. J. Phys. A: Math. Theor., 41:095003, 2008. 

[7] G. Hystad. Periodic Ising Correlations. PhD thesis, University of Arizona, Tucson, AZ, 2009. 

[8] B. Kaufman. Crystal Statistics. II. Partition Function Evaluated by Spinor Analysis. Phys. Rev., 76:1232- 
1243, 1949. 

[9] B. Kaufman and L. Onsager. Crystal Statistics. III. Short-Range Order in a Binary Ising Lattice. Phys. 
Rev., 76:1244-1252, 1949. 

[10] O. O. Lisovyy. Tau Functions for the Dirac Operator on the Cylinder. Communications in Math. Phys., 
255:61-95, 2005. 

[11] J. H. McClellan and T. W. Parks. Eigenvalue and Eigenvector Decomposition of the Discrete Fourier 

Transform. IEEE Trans. Audio and Electroacoustics, 20(l):66-74, 1972. 
[12] L. Onsager. Crystal Statistics. I. A Two-Dimensional Model with an Order-Disorder Transition. Phys. 

Rev. (2), 65:117-149, 1944. 
[13] J. Palmer. Planar Ising Correlations. Progress in Math. Phys., Volume 49, Birkhauser, 2006. 
[14] J. Palmer and G. Hystad. Spin Matrix for the Scaled Periodic Ising Model. To Appear in J. Math. Phys.. 

Preprint at arXiv:1008.0352v2 [nlin.SI]. 
[15] B. Simon. Representations of Finite and Compact Groups. Graduate Studies in Mathematics, Volume 10, 

AMS, 1996. 

[16] M. Spivak. Calculus on Manifolds. W. A. Benjamin Inc., New York, New York, 1965. 

[17] E. T. Whittaker and G. N. Watson. A course of Modern analysis. An Introduction to the General Theory of 

Infinite Processes and of Analytic Functions; with an Account of the Principal Transcendental functions. 

Fourth edition. Reprinted. Cambridge University Press, New York, 1962. 



36 



